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Abstract 



We study the semi-classical trace formula at a critical energy level for a /i-pseudo- 
differential operator on M" whose principal symbol has a totally degenerate critical 
point for that energy. We compute the contribution to the trace formula of isolated 
non-extremum critical points under a condition of "real principal type". The new 
contribution to the trace formula is valid for all time in a compact subset of M but 
the result is modest since we have restrictions on the dimension. 

Key words: Semi-classical analysis; Trace formula; Degenerate oscillatory 
integrals. 



1 Introduction. 

The semi-classical trace formula for a self-adjoint /i-pseudo-differential op- 
erator Ph, or more generally /i-admissible (see [19]), studies the asymptotic 
behavior, as h tends to 0, of the spectral function : 



where the Xj{h) are the eigenvalues of Ph- Here we suppose that the spectrum is 
discrete in [E — e,E-\-e], some sufficient conditions for this will be given below. 
If Pq is the principal symbol of Ph we recall that an energy E is regular when 
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Vpo{x, 7^ on the energy surface = {{x, ^) e T*W^ / po{x, ^) = E} and 
critical when it is not regular. 

It is well known that the asymptotics of (1), as /i tends to 0, is closely related 
to the closed trajectories of the Hamiltonian flow of p on the surface E^, i.e. : 

lim 7(E, h, cp) ^ {{t, e K X Ee / ^t{x, = {x, e)}, 

ft— »u 

where : 

and — d^po-d^ — d^Po-d^ is the Hamiltonian vector field attached to Pq. 

When ii^ is a regular energy a non-exhaustive list of references concerning this 
subject is Gutzwiller [11], Balian and Bloch [1] for the physic literature and 
for a mathematical point of view Helffer and Robert [12], Brummelhuis and 
Uribe [3], Paul and Uribe [17], and more recently Combescure et al. [8], Petkov 
and Popov [18], Charbonnel and Popov [7]. 

When one drops the assumption that E he a. regular value, the behavior of (1) 
will depend on the nature of the singularities of p on Eg which can be com- 
phcated. The semi-classical trace formula for a non-degenerate critical energy, 
that is such that the critical-set €{po) = {{x,^) e T*M" / dpo{x,^) = 0} is a 
compact manifold with a Hessian d'^po transversely non-degenerate along 
this manifold has been studied first by Brummelhuis et al. in [2] . They studied 
this question for quite general operators but for some "small times", that is 
they have supposed that the support of (f is contained in such a small neigh- 
borhood of the origin that the only period of the linearized flow in supp(<^) is 
0. Later, Khuat-Duy in [15,16] has obtained the contributions of the non-zero 
periods of the linearized flow with the assumption that supp(<^) is compact, 
but for Schrodinger operators with symbol ^^ + V{x) and a non-degenerate po- 
tential V. Our contribution to this subject was to compute the contributions 
of the non-zero periods of the linearized flow for some more general operators, 
always with (p of compact support and under some geometrical assumptions 
on the flow (see [4] or [5]). Finally, in [6] we have obtained the contributions to 
the semi-classical trace formula of totally degenerate extremum and the objec- 
tive of this article is to obtain a generalization when one drop the extremum 
condition. 

Basically, the asymptotics of (1) can be expressed in terms of oscillatory inte- 
grals whose phases are related to the flow of po on E^;. When (xq, ^o) is a critical 
point of po, it is well known that the relation Ker{dx,^^tixo, Co) ~ Id) ^ {0} 
leads to the study of degenerate oscillatory integrals. Here we examine the 
case of a totally degenerate energy, that is such that the Hessian matrix at 
our critical point is zero. Hence, the linearized flow for such a critical point 
satisfies dx,^^t{xo, Co) = M, for alH e R and the oscillatory integrals we have 
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to consider are totally degenerate. 

The core of the proof lies in establishing suitable normal forms for our phase 
functions and in a generalization of the stationary phase formula for these 
normal forms. The construction proposed for the normal forms is independent 
of the dimension but the asymptotic expansion of the related oscillatory in- 
tegrals depends on the dimension and on the order of the singularity at the 
critical point. This explains why the main result is stated with a restriction 
on the dimension. 



2 Hypotheses and main result. 

Let P/j = Op'^{p{x, ^, h))he a. /i-pseudodifferential operator, obtained by Weyl 
quantization, in the class of /i-admissible operators with symbol p{x, ^, h) ~ 
J2h^Pj{x:0- This means that there exists sequences pj e E^(T'*M") and 
RN{h) such that : 

Ph=Yl h'pjix, hD,) + ViV e N, 

j<N 

where RN{h) is a bounded family of operators on L'^iM."-), for h < ho, and : 

E^(r*R") = {a : T*W C, sup \d^d^a{x,C)\ < Ca,(3m{x,C), Vq;,/3 e N"}, 

where m is a tempered weight on T*]R". For a detailed exposition on h- 
admissible operators we refer to the book of Robert [19]. In particular, po{x, ^) 
is the principal symbol of Ph and pi{x, ^) the sub-principal symbol. Let be $t = 
exp{tHp^) : T*W T*W, the Hamiltonian flow of Hp^ = d^poA - d^po-d^. 

We study here the asymptotics of the spectral function : 

7(^0,/.)= E (2) 

Xj{h)€[Ec-£,Ec+e] 

under the hypotheses (Hi) to {H4) given below. 

{Hi) The symbol of is real and there exists £0 > such that the set 
Po^{[Ec - Eq, Ec + £0]) is compact in T*M". 

Remcirk 1 By Theorem 3.13 of [19] the spectrum a{Ph) fl [Ec — e,Ec + e] is 
discrete and consists in a sequence Xi{h) < \2{h) < ... < Xj{h) of eigenvalues 
of finite multiplicities, if e and h are small enough. 

To simplify notations we write z — {x, ^) for any point of the phase space. 
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(Hi) On — Po^{{Ec}), Po has a unique critical point zq — (xo,^o) O'^d 
near zq : 

AT 

Po{z) = + EPi(^) + Oi\\{z - Zo)\r'), k > 2, 

j=k 

where the functions pj are homogeneous of degree j in z — Zq. 

We have e C^{R). 

Since we are interested in the contribution of the fixed point Zo, to understand 
the new phenomenon it is suffices to study : 

1 r tE i 
^,,{E,, h) = — Tr j e'^0{t)ij^{x, hD,)exp{--tPhMPh)dt. (3) 

Here is a function of locahzation near the critical energy surface Eg^ and 
'0 G C^(T*]R") has an appropriate support near Zq. Rigorous justifications arc 
given in section 3 for the introduction of Q{Ph) and in section 4 for il)'^{x, hD^). 

In [6] it was proven that : 

Theorem 2 Under hypotheses (Hi) to (H^) and if Zq is a local extremum of 
the principal symbol po we have 

N 

7.o(^c, h) ~ /iX-(E A + C»(/i^)), ash^O, 

j=0 

where the ^ are some distributions and the leading coefficient is given by : 

1 / 2n-fc \ I I- 2n 

Ao,fe(<^) = ^(¥'(t + Pi(^o)),V j \Pkm--de, (4) 

^ ^ S2n-1 

with tzQ = max(t, 0) if Zq is a minimum and t^g — max(— 0) for a maximum. 

To obtain a generahzation of Theorem 2, when the critical point Zq is not a 
local extremum, we consider the classical hypothesis : 

{H4) We have Vpk ^ on the set Cp, ^ {9 e S^""^ / pk{0) = 0}. 

Remark 3 We would like to emphasize that, contrary to the case of a local 
extremum, the critical point zq is not necessarily isolated on Se^- This imposes 
to study the classical dynamic in a micro-local neighborhood of zq- Moreover, 
by homogeneity {H4) implies that Vpfc 7^ on the cone {{x, 7^ / Pfc(^) — 
0}. This allows to define the Liouville measure on the set Cp^. 

Then, the new contribution to the trace formula is given by : 

Theorem 4 Under hypotheses {Hi) to {H4) and if k > 2n we have : 
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where the leading coefficient is given by : 

^ ' ^ — ^ \ 1 f 2n 



Ao,fe(<^) = ^ 



. , , (27r)- 

V S2„-in{p^>o} 



+ (^\tfI-\ip{t+p,{zo)))j^ J ip.ier'^dey (5) 

If k — 2n a similar result is : 

^,,{E,, h) ~ h^-^\og{h)k{ip) + as /i ^ 0, 

with leading coefficient : 



where LVol is the Liouville measure attached to pk restricted to the sphere 

LVol(CpJ= J dL,,{e), 



(6) 



with dLp^{6) A depkiO) = dO on Cp^. 

Remcirk 5 In Eq. (5) and under the assumption on the dimension both terms 
are well defined since for the first term 1^1"^"^ is locally integrable and for the 
second term the integral is convergent. 



3 Oscillatory representation. 



Let be e 5(M) with (p e C^(M), we recall that : 

with Po ^(/eo) compact in T*M". By Proposition 3.13 of [19] the spectrum of 
Ph is discrete in 1^ for /i > small enough and e < Sq. Now, we localize near 
the critical energy E^ with a cut-off function G C^(\Ec — e,Ec + s[), such 
that = 1 near Ec and < © < 1 on M. The associated decomposition is : 

^{E,, h)^^^{E,,h)+^2{Ec, h), 

with : 

J^{E^,h)^ y: {1 - emh)M ^^^^\- \ (7) 
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MEch)^ E e{Xj{h))cp C^\ (8) 

The asymptotic behavior of ji{Ec, h) is classical and is given by : 
Lemma 6 -ii{E^, h) = e)(/i°°), as /i ^ 0. 
For a proof see e.g. [6] . 

Consequently, for the study of 7(-E'c, h) modulo 0{h°°), we have only to con- 
sider the quantity 72 (£'0 h). By inversion of the Fourier transform we obtain 
the identity : 

Since the trace of the left hand-side is exactly 72 (-^c, h), we have : 

72(£e, h) = ^Tr J e''^mexp{-^tPH)e{P,,)dt, (9) 

M 

and with Lemma 6 this gives : 

^{E,,h) = i-Tr I e''tcp{t)eM-ltPhMPh)dt + 0{h°°). 

Let be Uh{t) = exp(— ^P/^), the evolution operator. For each integer N we 
can approximate Uh{t)Q{Ph), modulo 0{h^), by a Fourier integral-operator, 
or FIO, depending on a parameter h. Let A be the Lagrangian manifold asso- 
ciated to the flow of po, i-e. : 

A = {(t, r, X, e, I/, ri) e r*M x T*W x T*R" : r = 0, (a;, = r;)}- 



Theorem 7 T/ie operator Uh{t)Q{Ph) is h-FIO associated to A, there exist 
UQ^il{t) with integral kernel in /(R^"'^-'-, A) and R^j^\t) bounded, with a L^- 
norm uniformly bounded for < h < 1 and t in a compact subset of W, such 
that UnmiPn) = C/g2(^) + h''B!i'\t). 

We refer to Duistermaat [9] for a proof of this theorem. 

Remark 8 By a theorem of Helffer and Robert, see e.g. [19], Theorem 3.11 
and Remark 3.14, Q{Ph) is an /i-admissible operator with a symbol supported 
in Pq ^(-^e)- This allows us to consider only oscillatory- integrals with compact 
support. 
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For the control of the remainder, associated to Rf^ '{t), we use : 

Corollary 9 Let be Qi e C^(M) such that Qi = 1 onsupp(6) ano? supp(Oi) C 
le, thenyN e N : 

TT{e{P,M^^^)) = ^Tr J 0{t)ei'^^U^';i{t)Q,{P,)dt + O{h''). 
For a proof, see e.g. [6]. 

If {xo,^o) G A and ii (p = (p{x, 6) e C^{R^ X M^) parameterizes A m a 
sufficiently small neighborhood U of (xo,^o) then for each Uh G /(M'^,A) and 
X G C^{T*'E.^'), supp(x) C U, there exists a sequence of amplitudes aj — 
aj{x, 9) e C^{R'' X M^) such that for all L e N : 

X'"{x,hD,)u,= hU{ajei^) + 0{h'^). (10) 

-d<j<L 

We will use this remark with the following result of Hormander (see [13], tome 
4, proposition 25.3.3). Let be {T,T,Xo,^o,yo, —rjo) e Agow, % 7^ 0, then near 
this point there exists, after perhaps a change of local coordinates in y near 
yo, a function S{t, x, rj) such that : 

Ht^x^y^v) ^ S{t,x,rj) - {y,ri) , (11) 

parameterizes Aflow In particular this implies that : 

{{t, dtS{t, X, rj),x, d^S{t, x, rj), dr,S{t, x, rj), -rj)} C Aflow, 

and that the function 5* is a generating function of the flow, i.e. : 

Mdr,S{t, X, v),v)- (x, d^Sit, X, rj)). (12) 

Moreover, S satisfies the Hamilton- Jacobi equation : 

dtS{t,x,r]) + po{x,d^S{t,x,rj)) = 0, 
S{Q,x,0 = {x,0. 

Now, we apply this result with (xcCo) = {yo,Vo), our unique fixed point of 
the flow on the energy surface Eg^. 

Remark 10 If = we can replace the operator by ei<-'«^)P^e-i<-'«^) 
with ^1 7^ 0. This will not change the spectrum since this new operator has 
the symbol p{x, ^ — ^i, h) and the critical point is now {xq, Ci) with ^ 0. 
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Consequently, the localized trace 72(-E'c, h), defined in Eq. (9), can be written 
for all e N and modulo 0{h^) as : 

To obtain the right power —d of h occurring here we apply results of Duister- 
maat [9] (following here Hormander for the FIO, see [14] tome 4, for exam- 
ple) concerning the order. An /i-pseudo-differential operator obtained by Weyl 
quantization : 

is of order w.r.t. 1/h. Now since the order of Uh{t)Q{Ph) is — |, we find that 



iP^{x,hD,)Uh{t)e{Ph) - E(27r/i)-'^+^ y a,(t,x,y,7;)et(5(*'-''')-<^'''))dy. 

j<N 

(14) 

Multiplying Eq. (14) by (p{t)eh*^'' and passing to the trace we find Eq. (13) 
with d — n and where we write again aj{t, x, rj) for aj{t, x, x, rj). 

To each element Uh of /(M'^, A) we can associate a principal symbol et'^(jprinc(wfe) , 
where 5" is a function on A such that ^dx — dS on A. In fact, if — I{aei'^) 
then we have S — — (poi~^ and aprmc{uh) is a section of |A|2 (g)M(A), where 
M(A) is the Maslov vector-bundle of A and \A\^ the bundle of half- densities 
on A. The half-density of the propagator Uh{t) can be easily expressed in the 
global coordinates {t,y,rj) on Agow If Pi is the sub-principal symbol of Ph, 
then this half-density is given by : 

t 

cxp(« J pi{<^,{y,-r]))ds)\dtdydi]\^. (15) 



This expression is related to the resolution of the first transport equation for 
the propagator, for a proof we refer to Duistermaat and Hormander [10]. 



4 Classical dynamic near the equilibrium. 



A critical points of the phase function of (13) satisfies the equations : 



E,= -dtS{t,x,0, 
x^d^S{t,x,0, 
^^d,S{t,x,0, 



$,(x,0 = (^,0, 
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where the right hand side defines a closed trajectory of the flow inside S^^. 
Since we are interested in the contribution of the critical point, we choose a 
function ip e C^(T*M"'), with ip = \ near zq, hence : 

l2{E,,h) = ^Tr j e''^mi^^(^,hD,>M-ltPhMPh)dt 

+^Tr J e''f'ipit)il-^lj^{x,hD,))exp{-^tPh)eiPh)dt. 

Under the additional hypothesis of having a clean flow, the asymptotics of the 
second term is given by the semi-classical trace formula on a regular level. We 
also observe that the contribution of the first term is micro-local. Hence this 
allows to introduce local coordinates near zq. To separate the contribution of 
Zq from other closed trajectories we use the following result on the classical 
dynamic. 

Lemma 11 For allT > there exists a neighborhood Ut of the critical point 
such that ^t{z) 7^ z for all z e ^7t\{-2o} and for all t e]-T, 0[U]0, T[. 

Proof. Since Zq is a degenerate critical point we have dHp^{zo) — 0. Hence, for 
all £ > we can find a neighborhood U of zq such that : 

\\dHp^(z)\\ <s,yze U. 

By a theorem of Yorke [21] we obtain that any closed trajectory in U has a 
period Tq that satisfies Tq > 27re~^. Thus for any T we can choose £t, and 
then Ut, such that To > T. ■ 

With supp((^) compact we can choose iIj such that Lemma 11 holds on supp('0) 
for all t G supp((^). Hence, on the support of ijj there is two contributions : 

1) Points (t,x,0 = (O,x,0 for G S^,. 

2) Points (t,a;, ^) = {t,zo) for t e supp(<^). 

The first contribution is non-singular for (x,^) ^ zq and can be treated, again, 
by the regular trace formula. Now, we restrict our attention to the second 
contribution and, since Zq is totally degenerate, we obtain : 

d^tizo) = exp(O) = Id, Vi. (16) 

The next homogeneous components of the flow are given, with {H2), by : 

ci^$t(^o)=0, Vi, Vje{2,..,A;-2}. (17) 

To obtain the next non-zero term of the Taylor expansion of the flow, we will 
use the following technical result : 

Lemma 12 Let be zq an equilibrium of the C°° vector field X and $t the flow 
of X. Then for all m e N*, there exists a polynomial map P^, vector valued 
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and of degree at most m, such that : 

t 

d^Mzo)(z'^) ^ d^Zo) I d^^s{zo)Pm{d^s{zo){z),...,d'^-^^s(zo)(z^-'))ds. 



For a proof we refer to [4] or [5] . 

Since d^ti^o) = Id, for all t, the first non-zero term of the Taylor expansion 
of the flow is given by : 

t 

d'-'M^o){z''~') = / d'-'Hp,{zo){z''-')ds = td'-'H,^{zo){z'-^). (18) 



Where the last identity is obtained using that d^~^Hpf^{zo) — d^~^Hp^{zQ). 
Moreover, with ci^p(^o) — 0, for the next term Lemma 12 gives : 

t 

d'=$,(zo)(^') = / d^H,,{z,){z'')ds^td^H,^^^{z,){z''). 



Remeirk 13 For higher order derivatives d^^tiz^), with j > k, there is two 
different kind of terms, namely : 

t 

I d^HM{z^)ds = td^HM(z^) = 0(\\z\\^), 



and terms involving powers of t, for example we have : 
t 

ld^+'-'H,,{zo){z^^'-\d'-'%{zo){z'-')) 



= *^d^^'-'H,,{zo){z^^'-\ d'-'H,,{zo){z'-')). 

This term is simultaneously 0{t^) and near (0,2;o). A similar result 

holds for other terms, which are 0{t^) for d > 2, by an easy recurrence. 

Lemma 14 Near Zq, here supposed to be to simplify, we have : 

S{t, X, - {x, + tE, = -t{Pk{x, + Rk+i{x, + tGk+i{t, X, 0), (19) 

where R,+,{x,0 = O(||(a;,0ll'+') and Gk+i{t,x,0 = O(||(a;,0ll'+'), urn- 
formly with respect to t. 
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Proof. By Taylor we obtain : 

We search our local generating function as : 

N 

s{t, X, = -tE, + (x, + E s{t, X, + o{\\ {x, r+'), 

where the S'j are time dependant and homogeneous of degree j w.r.t. {x,^). 
With the implicit relation a;, ^) = (a;, a;, ^)), we have : 

S{t, X, = -tE, + {x, + Skit, X, + 0{\\ {x, 1 1^=+^), 

where Sk is homogeneous of degree k w.r.t. (x,^). If J is the matrix of the 
usual symplectic form, comparing terms of same degree gives : 

By homogeneity and with Eq. (18) we obtain : 

Sk{t,x,0 = l^{{x,0,tJd'-'H,^{x,0'-') = -tPk{x,0, 
Sit,x,0 = -tE, + {x,0-tPk{x,0 + Oi\\{x,0\\'^'). 

As concern the remainder, we first observe that S{0,x,^) = (x,^). Hence, we 
can write : 

S{t,x,0-{x,0-tFit,^,0, 

with F smooth in a neighborhood of (x,^) = 0. Now, the Hamilton- Jacobi 
equation imposes that F{0,x,^) = —po{x,^) and we obtain : 

Rk+i(x,O^Po(x,0 - E,-Pk{x,0 = O(||(x,0ll'+')- 
Finally, the time dependant remainder can be written : 

S{t, X, - S{0, X, - tdtS{0, X, = Oit^), 

since by construction this term is of order 0(||(a;, (^)||'^"'"^) we get the desired 
result when i is in a compact subset of R. ■ 
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5 Normal forms of the phase function. 



Since the contribution we study is local, we can work with some coordinates 
and identify locally T*W^ with R^" near the critical point. We define : 

^{t, z) = ^{t, X, = S{t, X, - {x, + tE,, z = {x, e M^n. (20) 

Lemma 15 If Ph satisfies conditions {H2) and {H4) then, in a neighborhood 
of {t, z) = (0, ^0); there exists local coordinates x such that : 

"^{t, z) ~ — XoXi) ^'^ (^J-^ directions where Pk{0) > 0, 
^(t, z) ~ +X0X1) '^11 directions where Pk{0) < 0, 
^{t, z) ~ X0X1X2, near any point where Pk{0) = 0. 

Proof. We can here assume that zq is the origin and we use polar coordinates 
z = {r,6), ^ e §^'*~^(R). With Lemma 14, near the critical point we have : 

^{t, z) ~ -tr'iPkie) + rRk+i{e) + tGk+i{t, rO)), 

where pk{d) is the restriction of on S^""^. 
If pfc(^o) 7^ we define new coordinates : 

(XO, X2, X2n)(t, r, e) = (t, 01, 92n-l), 

Xi(t,r,^) = r\Pk{9)+rRk+i{e)+tGk+i{t,re)\-K 

In these coordinates the phase becomes — XoXi if Pfc(^o) is positive (resp. XoXi 
for a negative value). Near 9q, we have : 

|^(i,o,^) = M^)lMo, yt, 

hence, the corresponding Jacobian satisfies 17x1(^,0,6*) = |pfe(6')|^ 7^ 0. 

Now, let ^0 be such that pk{9o) = 0. Up to a permutation, we can suppose 
that de^pkiOo) ^ 0. We choose here the new coordinates : 

(X0,Xl,X3, ■■■,X2r^{t,r,Q) = (t,r,6'2, ...,6'2n-l), 

X2(t, r, 9) = Pk{9) + rRk+i{0) + tGk+i{t: r9). 
Since we have | Jxl(^) 0, 6*0) = \de^pk{9o)\ 7^ 0, lemma follows. ■ 
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In order to use these normal forms we introduce an adapted partition of unity 
on S^"~^. We choose functions ^j{0), with compact supports, such that : 

{9 e §2"-^ / Pk{9) = 0} c Usupp(^), 

j 

so that normal forms of Lemma 15 exist inside supp(f2j). Since is com- 

pact this set of functions can be chosen finite and we obtain a partition of 

unity with = 1 — The support of Qq might be not connected and 

we define Qq, with pk{0) > on supp(no similarly we define where 

pk < 0, so that Qq = ^0 + ^0 • 

If we accordingly split up our oscillatory-integral we obtain : 

7+(A)= J ei''^''''''^'^n+{0)a{t,r0y-^dtdrd0 

RxM+xS2n-l 

= J e-y^^°>^"A+{xo,Xi)dXodxi, 

MxM+ 

for the directions where pk{d) > and : 

7_(A)= J ei^^''''^^noie)a{t,r0y-^dtdrd0 

RxM+xS2"-i 

= j eT^^°^'Ao{xo:Xi)dXodXu 

MxM+ 

for the directions where pk{0) < 0. Similarly, the contribution of the neigh- 
borhood of the set {pk{d) — 0} is given by : 

/,(A)= J ei^'''-''''''^nj{0)a{t,r0y-^dtdrd0 
= / e^^°^f^Mj(xo,Xi,X2)c?Xoc?Xic?X2. 

RxR+xR 

The associated new amplitudes are respectively given by : 

^^(Xo,Xi) = / X*i^^iO)ait,r0y'~'\Jx\)dX2...dx2n, (21) 
^,(Xo,Xi,X2) = / x*{^j{OHt,r0y^~'\Jx\)dx3...dX2n. (22) 

Remark 16 Since Xi(t,r, ^) = r\pk{0) + rRk+i{0) + tGk+i{t,r9)\^ , our new 
amplitude satisfies ^o'(XoiXi) = ^{Xi"'~^)j near Xi = 0. A similar argument 
shows that A(xO)Xi)X2) = (^(Xi"^^)) near xi = 0. These facts will play a 
major role in Lemmas 17, 19 and 21 below. 
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We end this section with lemmas on asymptotics of oscillatory integrals with 
phases as in Lemma 15. 

Lemma 17 There exists a sequence {cj)j of distributions, whose support is 
contained in the set {xi — 0}, such that for all function a e C^(R+ xR) ; 

7 r °° +1 

J ( J e^^'^'^^'aixo, Xi)dXo)dxi - E >^~'~cj{a), (23) 
ffi J=o 

asymptotically for A — > oo, where : 

1 1 j+i-fc 

Cj = - — {J^{x_ )(xo) ® ^0 (Xi))> = max(-x, 0). 

We refer to [6] for a proof of this lemma. ■ 

Remark 18 A similar result holds for a phase —XoXi if we replace terms a;_ 
by x+ in Lemma 17. 

Lemma 19 If k > 2n we have : 

oo 

J ( I e'^^^^'^^^a{xo,Xi,X2)dxodX2)xf'-'dxi = A-^d(a) + 0(A-^log2(A)), 

R2 

where the leading coefficient is given by : 

d{a) = ^r(y) J \xoX2\'^exp{i^sign{xoX2))a{xo,0,X2)dxodX2- (24) 



Proof. We use the Bernstcin-Sato polynomial, see [20] for a detailed construc- 
tion. We use variables {t,r,v) instead of (xo,Xi)X2) and since Xi ^ foi" 
tv >0 we can write : 

:((tt;r'=)i-V"-^) = 6fc(z)(tOT^)-V2"-\ (25) 



dtdv dr^ 

k 

bk(z) = (1 - z)2 n(i - A;^ + 2n - 1). (26) 

3=1 

With the classical representation : 

oo 

J J e'^^'"^''a{t,r,v)dtdvdr 

{tv>0} 

= — y e^'?r(^)A-"( y y Hvr'^yait, r, v)dtdvdr)dz, 

7 {to>0} 
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where 7 —]c — ioo,c + ioo[ and Re(c) < k~^, we can compute the asymp- 
totic expansion with the residue method by pushing of the complex path of 
integration 7 to the right. Note that, when the phase is negative, we have : 



1 
2in 



J e-'^T{z)X-'{J J (\tv\r^)-'a(t,r,v)dtdvdr)dz. (27) 



7 {tv<0} 

With Eq. (25) and (26) the meromorphic extension is given by : 

00 

j J {tr''v)-'r^''-^a{t,r,v)dtdrdv 

{tv>0} 



J J {tvr'^r^r^--'^—^a{t,r,v)dtdvdr. 



{tv>0} 

Under our assumptions, the poles are ^ = 1 as a double pole and : 

j + 2n-l 



k 

2n 



Since k > 2n, the first pole is ^ ^ Z. The residue in this pole is : 

2n n 2n f \i_2n k■—^ d"^ , \ 1 1 1 

Ck,n\ ^exp{zn-)T{—) J {tv) 1^ r ■^-^a{t,r,v)dtdvdr, 

{tv>0} 



c,,„ = (-1)^ hm 1— -f^ = (-1)^ 



h{z) ' ' {k-2nyr{k) 
Now, with the following relations : 

d'^a , , , , ,^ 2n 



J {tvf'^^^{t,v)dtdv = [l- j {tv)~^a{t,v)dtdv, 

t,v>0 t,v>0 

1)'/ r'^-'^dr = {k-l)\a{Q), 



.id'' a 

r 

r>0 

and also using that : 



k 2n,2„ .^,T..2n, r(f) 



-(i--)^(^-i)!r(-) 



{k-2n)^r{k)^ k' ^ ^ k' k ' 

we find that the first residue is given by : 

exp(i7r^) , 2n ,2n, /• , 2,, , 
^\ fc^ A~~r — ) / {tv)--a{t,Q,v)dtdv. 
k k J 

{t,v>0} 
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The summation over the other quadrants gives the resuh with Eq. (27). Fi- 
nally, the remainder is of order 0{X fe~log^(A)) ii 2n + 1 — k and of order 
0{X ^) otherwise. ■ 

Remcirk 20 The preceding result holds again, in a weaker sense, if k is not 
a divisor of 2n, as can be shown by iterating the process of meromorphic 

extension above and using the fact that our amplitude is 0{xi"^^^) near Xi = 0. 
Now, if k divide 2n then terms A^''log''(A), with d = 1,2, can appear for the 
leading term, since we obtain poles of order 3 for the first non-zero residue. 

Example. 

If we consider Gaussian amplitude, we obtain easily : 



7(A, k,n)^ I e'^''''''e-^''+'''+'-\''''-^drdtdv ^ 2n f " ^'^ dr. 







The choice of A; = 5 and n — 2 leads to : 

/(A, 5, 2) - 27r^2)£iiilA-| + 0{\-'). 
25y/n 

This extra Gamma-factor comes from the identity : 

\tv\-le-''e-'''dtdv = r(^)2. 



Lemma 21 When k = 2n we have the particular result : 

oo 

R2 

with the leading coefficient : 

d(a) = -a(0,0,0). (28) 
n 



Proof. We use again the Bernstein-Sato polynomial method. We use variables 
{t,r,v) instead of (xo,Xi)X2) and we define : 

oo 

J+{X) = J J e'^''''''''a{t,r,vy-'dtdrdv, 

{ti.>0} 
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similarly, we define J- (A) on the set {tv < 0}. Then we can write : 

1 °° 
J+(A) = ^ / e''^r{z)X-' j j {tv)-'r-^'''a{t, r, vy-^dtdrdv, 

7 {tv>0} 

where 7 —]c — ioo, c + ioo[, c < {2n)~^. Similar computations as in the proof 
of Lemma (19) show that the associated Bernstein-Sato polynomial is : 

2n 

b2n{z) = (1 - f[{j - 2nz + 2n- 1), 

i=i 

under our assumptions the first pole is 2; = 1 and is of order 3. We define two 
holomorphic functions, near z = 1, via : 

b2n{Z) 

Hence, when using the residue method, the first terms of the asymptotic ex- 
pansion of J+(A) are given by the formula : 

V {tv>0} J 

and we have a similar formula for J_(A) if we use G~{z) and integration over 
{tv < 0}. For all holomorphic application / and all A > 0, we have : 

{f{z)X-^) = ^{z)X-^ _ 21og(A)^(^)A-^ + log2(A)/(^)A-^ (29) 



dz^ dz^ dz 

The term involving log^(A) is computed as in Lemma 19 and we have 

The distributional factor is here given by : 
r 

{tv>0} 

Hence the contribution is given by : 

-log^(A)A-^-a(0,0,0), 
n 

but when the phase is negative we have : 

limG';(z) = +^, 



I I ^"''-"Q^Q^<t,^^^)dtdrdv^2{2n-l)\a{W). 
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consequently, by summation, there is no term associated to log(A)^. 

Now, we compute the main term associated to the log( A) . 
For the same reason as previously all terms obtained by derivation of r(^) 
and {z — l)'V^2n(-2) will give a zero contribution by summation. This will 
not be true for terms obtained by derivation of the exponential. As concerns 
derivation of the meromorphic distributions, we obtain respectively : 



oo 



/ / \og{tvr'-y--'^-^a{t,r,v)dtdrdv, 



(2n)\ 

{tv>0} 



iff 

^ ' {to<0} 



By summation and using that log(|tf |r^") = logd^vj) + log(r^"') we easily 
obtain that the associated factor is given by : 



'2n 



\og{^v) a{t, 0, v)dtdv — J log(|fa|) ^^^ a{t,0,v)dtdv 

\tv>0} {tv<0} J 



A new splitting of the logarithms and integrations by parts show that the 
associated contribution vanish. 

It remains now to compute the term associated to the derivation of the expo- 
nential. An easy computation shows that the associated contribution is : 

^ ^ {tv>0} 

by integrations by parts, we finally obtain that the contribution is given by : 

{tv>0} 



\og{\)\-\i^f-^ J ^a{t,0,v)dtdv^-^log{X)X-'a{0,0,0). 
{tv: 

With Eq. (29) we have : 



J+(A) = ^ log(A)A-^a(0, 0, 0) + 0(A-^). 

A totally similar calculation gives : 

J-W = ^ log(A)A-^a(0, 0, 0) + 0{X-'), 

and the result holds by summation of J- (A) and J+(A). 
Example. 



18 



We use again an amplitude that is a product of Gaussian-functions, we have : 



,.xM2 

Hence, iov k — 2 and n = 1 we obtain : 



I{\ 2, 1) = i-(-27rVo(^) + (27r - log(l) + 7rHo(^)), 

where Jy(x), Yy{x) are respectively the standard Bessel functions of first and 
second kind. Also, H„(2;) is the standard Struve function defined by : 

zY{z)+zy\z) + {z'-v^)y{z)^-- 



7r(2v - 1)!!' 

From classical properties of these special functions, we obtain : 

/(A,2,l)~(ii2l(^)+0(A-). (30) 

Here, 7 is Euler's constant and is obtained by derivation of the T{z) factor in 
the formula that gives meromorphic extensions of our distributions. ■ 



6 Proof of the main results. 



Directions where pk{d) 7^ 0. 

Following step by step the proof of Lemma 4 of [6] we obtain that the first 
non-zero coefficient is obtained for Z = 2n — 1 (see Remark 16) and is given by 

11/ 2n-fc to ^\ \ 1 /■ 2n-k 

fe (2n-l)! V ^^^^ )®5o \^o^(Xo,Xi)) = ^y )(XoMo+(xo, O)dxo. 

Since by construction : 

4+(xo,o)= j aixo^mtmpkieT'-^de, (31) 

g2n-l 

we obtain that the local contribution, associated to supp(Qo is : 

)(xo),a(xo,o)) y n+{e)\p,{e)\--de. (32) 

S2n-1 

A similar computation gives the contribution of supp(no ), via : 

^((^(a;r^)(Xo),a(xo,0)) / %{e)\p,{e)\-'^de. (33) 

S2n-1 
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Now, since a{t,0) = 0{t) exp{itpi{zo)), cf. Eq. (15), the contributions of the 
directions where pk{d) 7^ are given, respectively, by : 

/+(A)~iA-^(|t|J^,^(t+pi(;.o))) / nUeMier'^de, (34) 

S2n-1 

for the set of directions where pk{0) > and by : 

7_(A)~iA-^(|t|!^,^(i+Pi(^o))) / n^{9)\p,{e)\-'^de, (35) 

§2n-l 

for the directions where pk{0) is negative. 



Microlocal contribution of the set pk{0) — 0. 
Case of k > 2n. 

Here we examine the contribution of terms : 



/i k 



2n. 2n f 2n ,7171 , \N 7 / n \ 7 7 / 7 221+1, 2/1 \ \ 

-y\-j^)h^ j \XoX2\ exp(^— sign(xoX2))A(Xo,0,X2)aXoaX2+C'(/i " log (/i)). 



but Lemma 19 shows that these are given by : 

1^.2n., 
k ^1 

Writing the delta-Dirac distribution as an oscillatory integral leads to : 

^r(y)/i^^ J 6*^^^1x0X2!"^ exp(i^sign(xoX2))^i(Xo,Xi,X2)c?XoC?XiC?X2C?^, 
with the amplitude : 

^i(Xo,Xi,X2) = j x*{^i{d)a{t,re)\Jx\)dx3-dx2n- 
Now, we return to the initial coordinates and with Xi = ^ this gives : 

In^)^'^^ I ^niXoX2){t,r,9)\-'^n,{9)eM^^^^^^^^ 
if we use that (xo, Xi, X2)(^, 0, 61) = (^,0,^^(6*)), we obtain : 

ir(^)/iif J \tpk{e)\-'-^exp{i'^sign{tpk{e)))ni{0)a{t,O)dtd0. 

MxS2n-l 

Hence, when pkid) is positive, we have : 

[ \t\-'-^exp{i^sign{t))<f{t)e''P'^'°^dt 

R 

l„,2n, /, , Tin. , , vrn, ^, , 

= T^i^)(\t\+ exp i— + \t\_ ' ex.p(-i—),ip(t)e'^P'^^°^ 

k k \ k k 
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Similarly, when pk{d) is negative we have : 



ir(^) / |t|-^exp(-z^sign(t))(^(t)e'*^'i(^°)(it 

rv rv J rZ 



k ^ k' 



.Tin 



t|+ exp(-i— ) + exp(i— ), (p{t)e 



k 



k 



Hence, by summation the contribution is : 



,7rn. 



+i sin ( ) ( 1 1 1 - ^ sign (t ) ,0{t) e''P' ^ 

rv 



If we use the classical relations : 



^(|a;|^)(0 = -2sin(— )r(A + l)|e| 



-A-l 



^(|a;|"sign(a;))(e) = 2^ cos(— )r(A + 
we obtain, after some manipulations, that the contribution is : 



\ §2n-l 



-de 



+ (\t\^ 



^-hign(t),^(t + p,(zo))) J Q,(^)|Pfe(^)r^sign(p,(^))d^ . 



We can split the integral with respect to d9 into two parts to finally obtain : 



1 2n 

—h fe 
k 



( 



^-1 \ ( 2n 



£Sl—1 



S2"-in{pfc>0} 



-"^de 



s2n-in{pfc<o} 



With Eq. (34) and (35), by summation on the partition of unity the main 
contribution to the trace formula is given by : 



k (27r)'* 



1^1+ ,^{t + Pi{zo)) 



iPkier'^de 



in{pfe>o} 



+ (\t\J-\^{t + p,{zo))) J |Pfe(^)|-^ 



d0 



S2"-in{pfc<0} 



(36) 



21 



And this proves the first statement of Theorem 4 for /c > 2n. ■ 

Case of k — 2n. 

Here the contribution is given by : 

I e^*(*'^'^)l],(^)a(t, rOy^-'dtdrdO = J e^>^«^>^^ A(xo, Xi, X2)dxodxidX2, 

RXR+X§2n-1 

and we recall that the amplitude is given by : 

A(xo,Xi,X2) = / x*mO)ait,rey''-'\Jx\)dX3-dx2n. 
With Ai — Xi^'^Ai, from Lemma 21 we know that we have : 

/ e^nXoxix2A.^^^^^^^^^)^ln-idxodxidx2 = — log(/i)A(0,0,0) + 0{h). 
J n 

Where, by construction, we have defined : 

With z = {zi, Z2, Z3), we write the amplitude as : 

A(0,0,0) = J e-'<''^^°'^^'^'^^Mxo,Xi,X2)dzdxodxidX2, 

and by inversion of our diffeomorphism we obtain : 

i.(0,0,0) = — ^ / e-'^''^^°'^''^^^^ni{e)a{t,re)drdtdedz. 
Using that (xo,Xi) = {'tj^)^ by integration w.r.t. {t,r,zi,Z2) we get : 

i.(0,0,0) = ^ [ e-^^3^^(°'°'^^Qi(^)a(0,0)d^dz3- 

By construction X2{t, 0, 6) = pk{6) and, since pk{6) is an admissible coordinate 
on supp(f)i), we can use the change of coordinates u = pk{d), this leads to : 

A(0,0,0) = ^a(0,0) j e-'''''ni{e)dLp^{e)dudz3, 
with dpk A dLpi^{9) — dO. If we introduce : 

Cp, = {ee §^"-1 / pu{9) = 0}, 

the sum over the partition of unity gives : 

Y: f nmdL,,{9) = LVol{CJ, 
i / 
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where LVol is the Liouville volume attached to pk{d) on S^" ^. Finally, since 
we have a{t, 0) — ip{t)exp{itpi{zo)) , the second statement of Theorem 4 holds. 
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